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We investigate the superfluid properties of the inner crust of neutron stars at finite temperature for 
different pairing functionals. We generalize the formalism adopted in article [l[ to include the effect 
of the temperature to calculate the specific heat of each given Wigner-Seitz cell. The calculations 
are done for two pairing forces, Gogny Dl and Vj ow -k, with finite range and a density dependent 
contact interaction. We compare in such a way the effect of the pairing strength and of the range 
on the thermal properties of the inner crust. 
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I. INTRODUCTION 

The study of thermal properties of the surface of neu- 
tron stars represents a possible way to provide strong 
constraints on its inner structure. In particular some au- 
thors have argued that the occurrence of rapid cooling 
processes in the interior of the neutron star and the ob- 
servation of the variations of the surface temperature can 
give a possible constraint on the equation of state (EoS) 
of the dense matter. It follows that it is important to 
obtain a precise model describing the cooling mechanism 
of the star. 

It has been shown 0, Q that the cooling time of the 
neutron star is mainly determined by the thermal re- 
sponse of the inner crust. This region of the star is formed 
by nuclear clusters surrounded by a sea of unbound neu- 
trons and ultra-relativistic electrons (see ref. @ and ref- 
erences therein). Thus the importance of performing ac- 
curate calculations of its specific heat, Cy. Within the 
literature, it exists already several calculations for the 
inner crust and for its thermal properties. In a recent ar- 
ticle, Baldo et al. [l3| has shown that pairing correlations 
play an important role for this region of the neutron star. 
Including such correlations the sequence of Wigner-Seitz 
(WS) cells is remarkably modified compared to the orig- 
inal one proposed by Negele and Vautherin [l4|- Several 
authors have done calculations in this system using non- 
relativistic functionals and with zero-range pairing inter- 
actions or finite range pairing interactions |ll], [T2J , 

but without performing a systematic comparison among 
the two. 

For this reason, we decided to continue the investiga- 
tion we started in our previous article [l[ (hereafter called 
article I) on the properties of different pairing functionals 
at finite temperature. In particular we decided to study 
three different pairing interactions: two with finite range 
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TABLE I: The WS cells representing different density regions 
of the inner crust and used in the present study. In the dif- 
ferent columns we have: the particle numbers Z, N, the WS 
cell radii -Rws, the background density p\ (obtained averag- 
ing the neutron gas density far away from the cluster) and 
its Fermi momentum k F . In the last column we also give k F , 
which is the Fermi momentum calculated using the average 
neutron density of the WS cell. 



{i.e. Viow-k and Gogny Dl), but with different strengths 
and one zero-range fitted to reproduce the properties of 
the finite range one in case of infinite systems. This al- 
low us to better comprehend the role of the range of the 
interaction and of the strength in our results. 

As done in article I, we adopt the Wigner-Seitz ap- 
proximation to study the different regions of the inner 
crust. For the validity of this approximation we address 
to the article of Chamel et al. [15( . Within each WS cell 
we perform fully self-consistent Hartrcc-Fock-Bogoliubov 
(HFB) calculations, in this way we can obtain the en- 
tropy and the specific heat of each cell as a function of 
the temperature for both protons and neutrons. 

The article is organized as follow in Secjlll we present 
the formalism of our calculations in both Wigner-Seitz 
cells and pure neutron matter (PNM), in Sec lIIII we dis- 
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Where q stands for neutrons (n) and protons (p); with 
/iF,q we indicate the chemical potential (also called Fermi 
energy) [TtJ ■ The latter plays the role of a Lagrange mul- 
tiplier to constraint on average the number of particles 
for each species separately. We used the standard nota- 
tion nlj for the spherical single-particle states with radial 
quantum number n, orbital angular momentum / and to- 
tal angular momentum j. ET^y and are the Bogoli- 
ubov amplitudes for the i-th quasiparticle of energy Efy . 
^n'nij 1S t ne thermal averaged single particle hamiltonian 
calculated using a Skyrme functional. In this case, as il- 
lustrated by Bonche et al. [H, HH , the structure of h q n , nl - 
is formally the same at zero temperature, what changes is 
the structure of local densities that now include an addi- 
tional Fermi-Dirac distribution. The explicit expression 
of such quantities can be found in ref. [7| . 



FIG. 1: (Colors online)We compare the value of the neutron 
pairing gap at the Fermi energy in symmetric and pure neu- 
tron matter calculated using the Vj ow -k (squares and circles 
respectively) and the DDDI interaction defined in Eq[5] (solid 
lines). See text for details. 

cuss the results we obtained concerning the densities, 
pairing gaps and specific heats and finally in Sec. IIV1 
we present our conclusions. 

II. CALCULATION DETAILS 
A. Inner crust of neutron stars 

We assume that the cluster structure of the inner crust 
is the one given by Negele and Vautherin [l4[. We sum- 
marize the results in Tab|TJ The cells are assumed to be 
spherical on the average and non-interacting among each 
others. 

The HFB equations are extended to include the effects 
of finite temperature and they take the form [l6| 



The HFB equations are solved in a spherical box us- 
ing the Dirichlet- Neumann mixed boundary conditions as 
done in the original Negele and Vautherin's article flij ]. 
As indicated in article I we have two choise for the bound- 
ary conditions : (i) even-parity wave functions vanish at 
R = Rws] (h) the first derivative of odd-parity wave 
functions vanishes at R = Rws ■ We call them Boundary 
Conditions Even (BCE), while Boundary Condition Odd 
(BCO) when we do the opposite. In the rest of the article 
we will always use the BCE condition, unless when explic- 
itly indicated the BCO ones. For the ph channel we adopt 
the Skyrme functional SLy4 [2(1 Hl| . The dependence of 
the results on the choice of the Skyrme functional has 
been already discussed in article I. For the pp channel 
we consider three different two-body pairing interactions: 

(i) a density-dependent contact interaction (DDDI) (22j; 

(ii) the finite-range Gogny Dl interaction 12311 : ( in) low- 
momentum realistic interactions (Viowfc) [24H27| . More 
details about the finite range pairing interactions have 
been already given in article I. 
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According to our previous results, we found that the 
simple DDDI, adopted in article I, is not adequate to 
treat on equal footing both neutron and proton superflu- 
idity. Following Margueron et al. [28|, H^j , we introduce 
a more complicated form for the contact interaction 
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where po{r) is the isoscalar density and p sa t = 0.16 fm 3 is the saturation density . The parameters /f ,/| are 
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To avoid the utraviolet divergency we introduce a 
sharp cut-off of 60 MeV on the quasiparticle energy [3(| ■ 
The strength of the interaction, V = -458.0 MeV fm 3 , 
is related to the cut off via the standard relation 1311 
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where a nn = —18.5 fm is the neutron- neutron scat- 
tering length. The parameters ryi = 0.71, 772 = 1-03, 
Qi = 0.51, a 2 = 0.51 are fixed to reproduce the results 
obtained with Vj w-k pairing interaction in both Symmet- 
ric Nuclear Matter (SNM) and Pure Neutron Matter us- 
ing the SLy4 mean field. In FigQ]we show the results of 
our fitting procedure. We notice that with our choice of 
parameters 77^,0^ the DDDI gives a satisfactory descrip- 
tion of the pairing gaps of the Viow-k interaction at the 
Fermi energy. It is important to recall that the DDDI in- 
teraction can reproduce the matrix elements of the V[ow-k 
interaction only at the Fermi energy, we refer to the dis- 
cussion in article 24j for more details. 

As done in article I, we drop the Coulomb term in the 
proton-proton pairing channel, this leads to an overesti- 
mate of about 10% of the proton pairing gaps. 

To have a better comparison with the existing liter- 
ature, in this work we include also the proton-electron 
interaction 133 that reads 
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This equation is valid under the assumption of uniform 
electron distribution. When presenting the results for 
the HFB pairing gaps, we show the Lowest-quasiparticle- 
energy Canonical State (LCS) pairing gaps [H, H3] • It 
is defined as the diagonal pairing matrix element cor- 
responding to the canonical single-particle state, whose 
quasi-particle energy 
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is the lowest. Here e q nl j stands for the diagonal matrix 
element of the single-particle field h q in canonical basis 
and the corresponding diagonal pairing-field matrix 
element. 



Before studying the effect of the presence of bound 
nucleons in the center of the WS cell, we analyze the 
supcrfluid properties of the pure neutron matter at finite 
temperature. For a given neutron density p n in PNM, 
the gap and number equations for a given temperature 
T, have to be solved simultaneously [331 ] 
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fi n is the neutron chemical potential, k B is the Boltz- 
man constant and E n {k) — (e n (k) — fJ, n ) 2 + A n (fc) 2 is 
the quasiparticle energy, while the single-particle energy 
e n {k) is given by the sum of the kinetic energy and the 
Hartree-Fock potential U^p 
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Where m* is the neutron effective mass. The tem- 
perature dependence is introduced through the term 

tanh ^ E 2t}gT ) • The number equation (cf. Eq. (|5])) pro- 
vides the relation between the density and the chemical 
potential [i n . In the limit of weak coupling, where A„ << 
£F,m the chemical potential can be approximated by the 
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entropy of the infinite system for unit volume is 
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where f n (k) = yL + exp j is the Fermi distribu- 
tion, and g = 2 is the degeneracy of the system. 
To test the accuracy of our code we compare in Fig|2] (a) 
the solution of the infinite system using EqfT] and using 
the HFB equations defined in Eq[TJ for a system without 
protons. The HFB equations, EqlU are solved in a box of 
35 fm with the BCE conditions. The test has been done 
adopting the V[ ow _k pairing interaction together with the 
Skyrme functional SLy4, that determines the neutron ef- 
fective mass, to*. We observe that the two methods are 
numerically equivalent. In FigJ2] (b) we show the en- 
tropy of the systems as a function of the neutron Fermi 
momentum, k F , calculated with the two methods. We 
notice again that the two techniques used to describe the 
PNM system are numerically equivalent. This remains 
true also for the others pairing interactions. 

It is possible now to define the specific heat per unit 
of volume as 
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FIG. 2: (Colors online) On the left, we show the neutron pairing gap Eq[7]in PNM for the Vi ow -k pairing interaction at different 
temperatures. On the right we show the density of entropy of the system as a function of the neutron Fermi momentum, hp. 
In both cases the solid lines correspond to the solutions of Eq[71 while the dots are the solution of Eq[T] See text for details. 



III. RESULTS 
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FIG. 3: (Colors online) Specific heat C v of neutrons, Eq llll 
calculated in PNM performing the derivative of the entropy 
calculated using Eq. 1101 See text for details. 
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The result is shown in Fig[3] for different values of the 
Fermi momentum as a function of the temperature. The 
specific heat increases as a function of the temperature, 
up to a given critical value, called T c , where it presents 
a sharp drop. Within the simple mean field approxima- 
tion, this temperature, T c , represents the passage from a 
superfluid regime to a non-superfluid one. This can be 
simply verified, by comparing the values given in FigOl 
from the one we can extract from Fig[5] (a) where the 
pairing gap goes to zero. 



A. Density 

In this section we discuss the results obtained by solv- 
ing the HFB equations, EqJTJ for the WS cells given in 
TabU As discussed already in article I, we do not show 
the results for the cell 982 Ge, since the calculations, also 
at finite temperature, do not converge toward a stable 
configuration. This problem of convcrgcncy could be due 
to the presence of a deformed minimum that is more en- 
ergetically favorable [34j . If we consider the occupation 
probabilities of some given Hartree-Fock level at zero 
temperature, they are just or 1. When we heat the 
system, this is no more true, since we introduce in the 
equations a Fermi-Dirac distribution that fragments the 
occupations of the levels (to some extent, similarly to the 
fragmentation introduced by a pairing interaction). As 
discussed recently in ref . [35| , this extra redistribution of 
the occupation of levels can favor the appearance of pair- 
ing correlations. We report in Appendix B, the example 
of 176 Sn; a double closed shell nucleus at T=0, that be- 
comes superfluid in a temperature range of [0.05, 0.63] 
MeV. 

In our calculations we work in the temperature inter- 
val [0-2] MeV, using a sampling step of fc^dT = 0.01 
MeV. As an example, we show in Fig|4] the variation 
of the density as a function of the temperature for the 
cell 1800 Sn calculated using the Mow-k pairing interaction. 
We observe that concerning the neutron, we do not have 
any remarkable variations going from fc^T = MeV to 
= 2 MeV. Since this variation of the occupations 
due to the temperature happens around the Fermi en- 
ergy, placed in the continuum, without any significant 
change on the neutron density. 

On the contrary, the proton density is more sensible to 
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FIG. 4: (Colors online) For the cell 1800 Sn, the neutron (solid 
line) and proton (dashed line) density profile at three different 
temperatures fc_gT = 0, 1, 2 MeV. The calculations are done 
using the Vi ow _k pairing interaction. In the inset we show the 
same result, but in semi- logarithmic scale. 

the effect of the finite temperature. At fcgT = MeV 
the proton density decreases exponentially, since the pro- 
tons occupy bound states, but with the increase of the 
temperature the proton can occupy scattering states. At 
ksT = 1 MeV the density of the proton gas is about 

p gas _ 10 -11 fm -3 and it readles p gas _ at 

ksT — 2 MeV. We can calculate the number of protons 
in the gas region using the formula 

N° as = [ Rn S ^r 2 dr [p T p {r) ~ p T p =\r)} , (12) 
JRo 

where Ro is the starting point of the integrations and 
is chosen to be R rj lOfm, pE{r) is the proton den- 
sity at finite temperature and pp" =0 (r) is the proton den- 
sity at zero temperature. Since we are interested only 
to the protons that occupy scattering states, in Eq|T2l 
we subtract the tail of the proton density at zero tem- 
perature. In such a way the result is not very sensitive 
to the choice of Ro. We thus obtain for k B T — lMcV 
the value of N° ut w 0.008, while for k B T = 2MeV we 
have N° ut sa 0.05 protons. A detailed analysis about 
the appearance of a proton gas for finite nuclei at finite 
temperature has been already presented by P. Bonchc 
et al. [TH, showing that the vapor is present because 
we impose our system to be confined in a box, otherwise 
the Coulomb repulsion should provoke the evaporation of 
this gas. Thus to have box-independent results, Bonche 
et collaborators suggested a vapor subtraction procedure. 

Since we work within the WS approximation, the pos- 
sibility of the presence of a proton gas is admitted, since 
the confinement is imposed by the presence of the other 
surrounding cells, although in such case the WS approxi- 
mation itself could be less justified. For the temperature 



regime in which we work, this proton gas is very weak 
and we will not apply any specific procedure , but we re- 
mark that this phenomenon should be better investigated 
in the future. 



B. Pairing gaps 

We discuss here the results we obtained from our HFB 
calculations at finite temperature, with special attention 
to the pairing gap. In FigJS]we show the neutron and pro- 
ton pairing gaps A'^g using the three different pairing 
interactions: Vi Q w-k (panel (a)-(b)), Gogny Dl (panel 
(c)-(d)) and the contact interaction (panel (e)-(f)). In 
all the figures of this article we present the results from 
180 Zr to 1800 Sn for all the pairing interactions plus the 
results of 1500 Zr for the Gogny Dl, since only with this 
interaction the results converge to a stable solution (see 
discussion in article I). 

We observe that as a common feature, the value of the 
LCS gap decreases as we increase the value of the temper- 
ature and it finally disappears when reaching a specific 
value of the temperature, called Tj?. A similar behavior 
has been already observed for the PNM system in Figj2] 
(a). Compared to the case of pure neutron matter, we 
observe that the A£ cs is not a smooth function of the 
temperature fcgT, but it presents some jumps. This is a 
consequence of Eq[6] that selects the pairing gap of the 
single-particle canonical state with the energy closest to 
the Fermi energy. As a result, it can happens that we do 
not always select the same level for different values of the 
temperature, since there could be a change of the single 
particle structure due to the temperature and the jumps 
correspond to the selection of a different state in EqJBl 
Another definition of pairing gap averaged on more levels 
around the Fermi energy would avoid this problem [i.e. 
see Eq |Bl[) . but it would have as inconvenient to mix the 
superfluid properties of the bound states, with the one of 
the external gas. For the case of protons, the two defini- 
tions of pairing gap give very similar results, since they 
are bound. 

A more detailed analysis on the properties of the gap 
LCS will be the subject of a forthcoming article [36j], but 
we can anticipate here that Eqj5] is a good filter in the 
case of WS cells to study the superfluid properties of 
the external gas, where we have usually the majority of 
matter in each cell. 

Comparing the panel (a) with the panel (e) of Fig. [5] we 
observe that the zero range force reproduces quite nicely 
the global trend and the absolute values of the neutron 
gaps as a function of the temperature, with a discrepancy 
of less than «10% on the values of the gaps obtained with 
the two forces. That is not the case for the protons, where 
there is almost a factor of 2 difference in the gaps, but 
both interactions predict a superfluid solution for protons 
(on the contrary of the DDDI adopted in article I). 

The origin of the difference comes from the fitting pro- 
tocol we adopted to determine our density dependent in- 
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FIG. 5: (Colors online) From the top to the bottom, neutron (left) and proton (right) LCS gaps, A. q LCS , as a function of the 
temperature and for the different cells given in Tab[I] In panel (a)-(b) we show the results for the Vi ow -k interaction, in the 
panels (c)-(d) the results for the Gogny Dl interaction and finally in panels (e)-(f) the results for the contact interaction DDDI. 
For the different panels, we adopt the same color code. See text for details. 



teraction (EqJ5]), since we used only infinite matter prop- 
erties (i.e. the value of the pairing gap in SNM and PNM) 
without taking into account any kind of shell effects (i.e. 
including some finite nuclei gaps into our % 2 minimiza- 



tion to determine the parameters 77^=1,2, cii=i, 2)- We ob- 
serve that the analytic formula obtained from BCS the- 
ory HH2I 
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FIG. 6: (Colors online) Critical temperature for neutrons 
feT" calculated using the three different pairing interactions 
discussed in the article. The solid lines represent the pure neu- 
tron matter result calculated using Eq[7] while the points rep- 
resent the WS calculations calculated using the BCE. On the 
same figure the crosses represents the calculations for three 
WS cells using Vi ow -k pairing interaction and BCO conditions. 
See text for details. 



k B T? « 0.56 • AI bT=0 , (13) 

used to derive the critical temperature from the value of 
the gap at zero temperature, A| bT=0 , is compatible with 
our findings. 

To understand the role of the bound nucleus in the 
center of the WS cell, in Fig|6]we compare the value of 
the critical density T™ at which the gap A r l cs for neu- 
tron disappears, in both PNM {i.e solving EqJ7|) and WS 
cell. We associate at each WS cell a neutron Fermi mo- 
mentum, hp, according to the value of the density of the 
external neutron gas, p b n , as given in TabJU We observe 
that similarly to the zero temperature case (see article 
I) , the presence of the cluster reduces the neutron super- 
fluidity for the high density cells. It can be seen from 
a smaller critical temperature for the WS cell compared 
to the PNM case. In Fig|6] it is not easy to clearly see 
such difference, but it will appear more clear when we 
will compare the specific heat of the two systems as done 
in Fig|8] and Fig|TUl we thus refer to the discussion in 
the next section. The Gogny Dl interaction presents the 
same qualitative behavior with the difference that usually 
the critical temperature, T c , appears at higher tempera- 
tures due to the stronger pairing strength compared to 
the Viow-k interaction (or the DDDI). This is what one 
can expect by simply using Eq ll3l for an estimation. We 
refer to Appendix A for a more detailed discussion. 

As a last remark, we observe that the DDDI interaction 
given in Eq(2] gives a very similar value of critical temper- 
atures for neutrons as the Vjow-k pairing interaction, as 
expected from the results of Fig[51 at most with a differ- 



ence of «30 KeV. Such small difference can not be clearly 
seen on this picture, but it becomes much more clear in 
FigJ7{a)(b), where we present the specific heat. In some 
cells and for few values of temperature T, our numerical 
code was not able to converge the HFB equations using 
the DDDI pairing interaction. The code oscillates among 
two possible solution that are relatively close in energy 
( AE sal MeV). We tried changing the damping param- 
eter [38| without any remarkable improvement. In this 
case we do not include these points in our calculations of 
the entropy and of the specific heat. 

We also tested the dependence of our results on the 
choice of the boundary conditions. As shown already in 
article I and also in ref.[39j], the BCE and BCO boundary 
conditions give similar results concerning gaps at &bT=0 
MeV. We checked that is true also at finite temperature. 
We performed calculations in 3 different cells: 200 Zr, 500 Zr 
and 1800 Sn using the Vi ow ^k paring interaction and BCO 
conditions. In Fig|6]wc report the result of such calcula- 
tions for the case of critical temperature. The differences 
among the critical temperatures obtained with the BCE 
and BCO conditions is at most «10 KeV in the cell 500 Zr, 
showing that the results are robust and do not depend 
on the particular choice of boundary conditions. This is 
in contrast with the results shown in ref . jT^| , the main 
difference among the two calculation is that Baldo and 
collaborators minimized a new set of WS cells using the 
different boundary conditions, while in this work we limit 
to the WS cells calculated by Negele and Vautherin [lj| ■ 
We also refer the reader to Appendix C for further discus- 
sions on the finite-size effects and boundary conditions. 



C. Specific heat 

Exactly as we did for the PNM system, we now ex- 
amine the entropy and the specific heat for the WS cell. 
The expressions for these quantities have been already 
given in EqllOl and Eq llll for the infinite system and it 
straightforward to adapt them for the calculations in the 
WS cells. 

In Fig[7] we show the specific heat of protons and neu- 
trons calculated using the Mow-k interaction. In each cell 
the proton specific heat is usually two orders of magni- 
tude smaller then the neutron specific heat. As already 
observed in FigJ3J the sudden variations in the Cy are as- 
sociated with a phase transition from superfluid to non- 
supcrfluid system. This can easily understood comparing 
the position of the peaks with the values of the critical 
temperature T£ given in FigO In the same figures in 
panels (a) (b) we also show the specific heat C v calcu- 
lated using the contact interaction given in Eq(2] We 
observe that the curves follow quite closely the one ob- 
tained using the V) w-k interaction, apart from a small 
difference for the high density cells. This is a direct con- 
sequence of the differences we observed in the Figsj5]and 
[5J since the DDDI interaction is able to mimic the super- 
fluid properties of the Mow-k only within a certain error. 
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FIG. 7: (Colors online) On the top row (a-b) the neutron specific heat C™ as a function of the temperature for the different 
cells of Tab U obtained with Vi ow -k pairing interaction, on the bottom (c-d) the specific heat for protons . In panel (a-b) 
we also show, as dashed-dotted line, the neutron specific heat C™ obtained using the DDDI pairing interaction given in Eq[2] 
See text for details. 



In FigJ3 (a), we observe that the low-density cells 
present two peaks. To better understand this phe- 
nomenon in Fig(51 we compare the neutron specific heat 
C v for the cell 180 Zr and the specific heat of PNM calcu- 
lated at the neutron Fermi momentum k F = 0.13 fm _1 . 
See Tabd 

The two systems show a similar behavior up to the 
critical temperature, feT™ « 0.1 MeV, where we have 
the first phase transition. For PNM the peak is located 
at hsT™ = 0.08 MeV and it corresponds to the phase 
transition from superfluid to non-superfluid system (as 
discussed in Fig(3]), then the C v monotonically increases. 
On the contrary for the 180 Zr cell we have a first peak 
at fcsT™- 1 = 0.09 MeV and a second one at k B T^ 2 = 
0.32 MeV. A similar result has been already discussed in 
ref.Q, adopting a simple contact force with no isovector 
dependence {i.e. putting /f = in Eq[5]). 

We can deduce that the first transition refers to the 
disappearance of the neutron superfluidity for the exter- 
nal neutron gas, fc^T 1 ™' 1 , while the bound neutrons in the 



cluster are still superfluid. The second peak at fcfiT"' 2 
corresponds to the transition of the bound neutrons from 
superfluid phase to non superfluid phase. To confirm this 
hypothesis, we analyze the local pairing field, A£ oc (i?), 
of neutrons for 180 Zr. We refer to article I for the exact 
definition of this quantity (see also ref. [4(J ) . 

In FigHwe show A£ oc (i?) for the WS cell 180 Zr at 
different temperatures calculated using the Vj ow -k pair- 
ing interaction. We clearly see that at the first critical 
temperature, /cbT™' 1 , the external neutron gas becomes 
non-superfluid, since the pairing field goes to zero, but it 
remains non-zero inside the nucleus without significative 
changes. At the second critical temperature, fc^r™' 2 , the 
superfluidity disappears also from the nucleus, and the 
entire cell becomes non-superfluid. We conclude that 
the value of the critical temperature of neutrons given 
in Fig|6] corresponds to the phase transition of the gas, 
since as we discussed previously in the text, we used as 
a probe to study the superfluid properties the gap LCS, 
A£ cs , given in EqjSJ As anticipated in previous section, 
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FIG. 8: (Colors online) The neutron specific heat for 180 Zr 
calculated with the Vj ow _k pairing interaction (solid line), and 
with the DDDI pairing interaction (dashed-dotted line). In 
the same figure we show the corresponding PNM calculations 
(dashed line) at k F = 0.12 fin -1 (see Tab0. 



FIG. 10: (Colors online) The neutron specific heat for the 
two cells 13 50-i800g n ^ golid y ine ^ and the corresponc iing PNM 

calculations (dashed line) at kp = 0.6,0.8 fm -1 (see Tab(I]) 
calculated using the V[ ow _k pairing interaction. In the same 
figure we show with the same color code, the results of WS 
cells calculated using the DDDI interaction (dashed-dotted 
line). See text for details. 
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FIG. 9: (Colors online)Local pairing field of neutrons, 
A2 oc (r), for tne WS cell 180 Zr at different temperatures 
calculated using the Viow-k pairing interaction. See text for 
details. 

this is a feature of the gap LCS valid for the WS cell 
system. 

According to FigJT] (a) we have 2 phase transitions in 
the cells: 180 Zr, 200 Zr and 250 Zr. As expected we find the 
same behavior, with the two phase transitions when using 
the DDDI interaction, confirming that this is a general 
property of the system. In Fig llOl we repeat the compar- 
ison already done in FigJSJ but for the cells 1350 Sn and 
1800 Sn. Again the DDDI and the Vi ow _k interactions give 
very similar results. We do not observe any major qual- 
itative difference when adopting the Gogny Dl pairing 
interaction, we thus refer the reader to Appendix A for 
a more detailed discussion. 



We now analyze the proton specific heat, C v , calcu- 
lated using the Mow-k pairing interaction. We observe 
that for the cells from 11 to 7 we have only one tran- 
sition occurring at temperature fc^X^ ks 0.8 MeV, see 
FigEl (c), that are usually at higher temperature than 
the corresponding phase transitions for neutrons. More 
interesting is the behavior of proton specific heat for 
the last three cells. For example in 1100 Sn we observe 
a first peak at /cbT)?' 1 = 0.62 MeV, corresponding to 
the passage supcrfluid/non-supcrfluid regime for protons 
(i.e. zero proton pairing energy). We encounter a second 
peak at Tf - 1 = 1.08 MeV this is the same value of the 
temperature of the neutron phase transition from the su- 
perfluid to non-superfluid phase k B T?' 2 = 1.08 MeV. We 
conclude that the modification of the neutron entropy af- 
fects the proton entropy trough the isovector interaction 
among neutron and protons at the mean field level. Since 
the number of neutrons is much bigger then the number 
of protons, we do not observe the opposite phenomenon 
(i.e. a change in the proton entropy affecting the neu- 
tron entropy). The results for neutrons obtained using 
the DDDI interaction arc in good agreement with the re- 
sults employing the Vjow-k pairing interaction as seen in 
FigiSlfTOl This is not the case for protons, where the two 
interaction gives quite different quantitative results. 

IV. CONCLUSIONS 

We have performed extensive calculations of Wigncr 
Seitz cell adopting three different pairing interactions: 
two pairing interactions with finite range, but with dif- 
ferent strengths, i.e. Vi ow -k and Gogny Dl, and a density 
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dependent contact one fitted to reproduce the infinite 
matter results given by the V[ow-k interaction. Since we 
used the same numerical code for all the forces, we have 
been able to make a strict comparison between them. 
Comparing with the results obtained in pure neutron 
matter, we investigated the role of the strength and of 
the range of the pairing interaction on the specific heat 
of the inner crust. As a result, we observed that the 
features of neutron superfluidity described by the V| ow _k 
interaction and our DDDI are essentially the same. This 
is not the case for the protons, where the two interac- 
tions are in agreement only at the qualitative level. We 
also investigated the effect of the strength of the pairing 
functional on the superfluid properties of the WS cells. 
In particular we confirm, as in article I, that this effect is 
more important for the high density cells, where pairing 
can also affect the configuration of the cell 0, Q~3| . We 
observe that, independently of the adopted pairing inter- 
action, the cells close to the outer crust shows two distinct 
phase transitions: the first one appears at low tempera- 
ture and it concerns the external neutron gas, while the 
second at higher temperature describe the transition of 
the nucleus from superfluid to non-supcrfluid phase. We 
think that the intermediate phase among the two critical 
temperatures, i.e. a, superfluid nucleus immersed in a sea 
of non-superfluid neutrons should be better investigated, 
since it could have important effects on the properties of 
the inner crust. In fact in this case we have a significant 
deviation among the specific heat of the WS cell and the 
infinite system at the same density. This investigation 
goes beyond the scope of the present article and we leave 
it for a future work. 

We conclude by observing, that our results are ob- 
tained within the mean field approximation. It has been 
observed by some authors |Tl| - t44| . that the inclusion of 
higher order diagrams, i.e particle- vibration coupling, 
can strongly affect the shell structure and the superfluid 
properties of the system. Such kind of calculations are 
still not available for WS cells, mainly because of the 
computational time requested. According to our findings, 
a possible alternative could be represented by perform- 
ing calculations in the infinite systems of both collective 
excitations (45l - l47j and pairing gaps [48[. Such results 
could be then used to fit an effective zero range pairing 
interactions. Finally the resulting contact force could be 
used to calculate the properties of the inner crust. This 
method although the underlying approximations, could 
help giving a first quantitative estimate of such higher 
order effects in WS cells. 
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Appendix A: Gogny Dl 

In FigQTJ we show the results of the specific heat for 
both protons and neutrons obtained using the Gogny Dl 
interaction. These results should be compared with the 
ones in Figf7]for the case of Mow k pairing interaction. As 
a general remark, we observe that compared to the Vi ow -k 
or DDDI results, we have the same qualitative results. In 
more detail, for the cells from 11 to 8 the specific heat is 
quite similar in the two calculations: in fact in this den- 
sity region, the two forces give similar values of neutron 
paring gap, see for example FiglHJ For higher density 
cells, the difference among the critical densities obtained 
with the Gogny Dl interaction and the Vi ow -k increases. 
In 500 Zr the two interactions give the same critical tern- 
peraturc k B T?_ Gogay m = /c B T"_ yiow k , while in 1800 Sn 
we have k B T?_ Go ^ D1 - k B T?_ v ^ k = 0.34 McV. This 
difference in the critical temperatures is simply due to 
a difference in the pairing strength and it can be under- 
stood using Eq |13l and comparing the values of the gaps 
at zero temperature for the two interactions in Fig. 5 of 
article I. 

It is interesting to observe that, even for this interac- 
tion, we have a two peak structure in the neutron spe- 
cific heat for the cells from 180 Zr to 500 Zr. Concerning 
the proton specific heat, we have a very interesting result 
for the cell 1500 Zr. In this case we have a clear two peak 
structure is shown in FigQT] (d) . As already discussed in 
the text, we observe that the position of the second peak, 
corresponds to the position of a peak in the neutron spe- 
cific heat in Fig[JTJ (b). Moreover, this is the only case 
where C v is only w j^Cy, when usually the ratio among 
the two is much smaller pa j^. We leave the investiga- 
tion of such cell for a future work, in fact according to 
the results of article I, it could be that 1500 Zr is stable 
at k B T = Mev and it becomes unstable (see discussion 
about Fig. 2 of article I) when the pairing gap vanishes 
around the critical temperature. 
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Appendix B: Pairing reentrance 

In this Appendix we shortly discuss the phenomenon of 
the pairing reentrance, recently discussed by Margueron 
et al. [35| . To have a better comparison with the exist- 
ing results already present in the literature, we perform 
calculations of 176 Sn using the SLy4 functional and the 
Mow-k pairing interaction. We define the pairing gap as 
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FIG. 11: (Colors online)We show the neutron specific heat C™ as a function of the temperature for the different cells of Tab|T] 
obtained with Gogny Dl pairing interaction, on the bottom the specific heat for protons. See text for details. 
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We observe from FiglEJ that at k B T = McV, the sys- 
tem is non superfiuid, then at fegT™ = 0.05 MeV we have 
a first phase transition: the nucleus becomes superfiuid 
in neutrons. The maximum value of the gap is reached 
around k B T w 0.45 McV and then at k B T? = 0.62 McV, 
the system return non-supcrfluid. This is quite different 
behavior compared to standard nuclei as see in ref. [49j or 
in our WS calculations shown in Fig. [SJ 
As discussed in the text, the Fermi-Dirac distribution, 
induced by the temperature, can provoke the rearrange- 
ment of the occupation of single particle levels around 
the Fermi energy. In this case, some particles are moved 
from the close shell to some resonant states in the contin- 
uum, thus allowing the possibility to turn on the pairing 
field in the last shell that is no more completely occupied 
as it was at zero temperature. We refer to Margueron 
et al. [35i | for a more detailed analysis of such results. 



As a last remark, wc underlined that for the case of an 
isolated nucleus, it would be important to project on the 
good particle number. Unfortunately, such method can 
not be used with most of the existing Skyrme functionals, 
we refer to ref. [50M5^ | for a more detailed discussion on 
this problem. We leave to a future work, a better analysis 
concerning the phenomenon of pairing rcentrance. 

Appendix C: Finite size effects 

We now shortly discuss how the discretization induced 
by the box affects our results. Since the main observable 
we have calculated in this article is the specific heat, we 
decided thus to compare the specific heat for the PNM 
case into different boxes of different size and compare 
the results with the one obtained by solving the HFB 
equations in the continuum (see discussion in Sec lIIB[) . 
We thus calculate here different WS cells without pro- 
tons. The value of the specific heat C v should not de- 
pend on the size of the box, so if wc observe differences 
among the two methods of doing the calculations, this 
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FIG. 12: (Colors online) Averaged neutron pairing gap 
Eg IB II for 176 Sn as a function of the temperature fcsT. 
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should originate from the discretization imposed by the 
box, similarly as it has been done in ref. [53| . The cal- 
culations have been done using the Skyrme SLy4 and 
Mow-k pairing interaction. In Fig |13l wc compare the spe- 
cific heat of neutrons for two types of calculations: (i) we 
solve the HFB equations given in Eq|7] for a given value 
neutron Fermi momentum, k F as given in TablU (ii) we 
solve the HFB equations in a box, EqfTJ for a fixed value 
of Rws without protons and using a number of neutron 
that reproduces the same value of k F in the given box, as 
described in TabU We observe that the two calculations 
are in good agreement for boxes bigger of ~ 25fm. When 
we use smaller boxes we observe the first difference due 
to discretization, it is the case of the cell Rws = 19-6 fm 
and kp = 1.08 fm -1 , where we have a small difference 
among the critical temperatures calculated with the two 
methods of about « 80 KcV. Very different is the case 
R ws = 14.4 fm and k F = 1.37 fm -1 . We observe that 
the behavior of the neutron specific heat is very different 
in the two cases, and we restore the agreement among the 
two calculations only in the limit of relatively high tem- 
peratures. Such calculations show clearly the limits of the 
WS approximations for high density cells. For such kind 
of calculations we also tested for the cell Rws = 27.6 
fm and k F = 0.81 fm -1 the dependence of our results on 
the different boundary conditions when we solve EqQ] in 
a box. In such case we observe that using the BCO or 
BCE does not affect the results on the specific heat and 
on the critical temperature. 



FIG. 13: (Colors online) Neutron specific heat, C v (in semi- 
logharitmic scale), as a function of the temperature for dif- 
ferent WS cells without protons. The calculations have been 
done (i) solving the HFB equations of EqQ (lines) and the 
(ii) HFB equations in a box (dots) , Eq[T] See text for details. 
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